Statistical properties of the truncated state with random coefficients 
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In this paper we introduce the truncated state with random coefficients (TSRC). As the coefficients 
of the TSRC have in principle no algorithm to produce them, our question is concerned about to 
what type of properties will characterize the TSRC. A general method to engineer TSRC in the 
running-wave domain is employed, which includes the errors due to the nonidealities of detectors 
and photocounts. 
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I. INTRODUCTION 



The generation of new states of the light field turned 
out to be an important topic in quantum optics in the 
last years, such as quantum teleportation 111, quantum 
computation yL quantum communication Q, quantum 
cryptography Q , quantum lithography decoherence 
of states [g, etc. The usefulness and relevance of quan- 
tum states, to give a few examples, cover the studyof i) 
quantum decoherence effects in mesoscopic fields 0; ii) 
entangled states and quantum correlations [8|; interfer- 
ence in the phase space ; collapses and revivals of the 
atomic inversion [lOj; engineering of (quantum states) 
reservoir |ll|: etc. Also, it is worth mentioning the im- 
portance of the statistical properties of one state to de- 
termine some relevant properties of another |l2j| as well 
as using specific quantum states as input to engineer a 
desired state [l3| . 

Here we introduce a truncated state having random 
coefficients (TSRC). This state is characterized by a con- 
stant phase relation between the relative Fock states, 
but with coefficients describing probability amplitudes 
obtained through some random number generator. Fea- 
tures of this state are studied by analyzing several of its 
statistical properties. We believe the TSRC can be a use- 
ful tool to test sequence of numbers whose generation is 
unknown from the observer. Indeed, preliminary results 
have corroborated this expectation [Tj| . 

This paper is organized as follows: in Section II we 
introduce the TSRC and Section III we analyze the be- 
haviour of some statistical properties of the TSRC when 
the Hilbert space is increased. In Section IV we show how 
to engineer TSRC in the running-wave domain, and Sec- 
tion V is devoted to the study of losses due to nonideali- 
ties of the detectors in the whole process of engineering. 
Finally, in Section VI we present our conclusions. 



Quantum state with random coefficients 



II. TRUNCATED STATES WITH RANDOM 
COEFFICIENTS (TSRC) 



We define the TSRC as 



N 



\TSRC) =AfJ2 r ne l 



(1) 



inO 



where J\f is a normalization constant and C n = r n e 
is a random complex coefficient (polar form), given by 
some random number generator (RNG). It is worthwhile 
to note that TSRC differs from (mixed) thermal states 
(MTS) and Pure States Having Thermal Photon Distri- 
bution (PSTD)H3, since these states have an algorithm 
for the coefficients C„/s. In addition, although TSRC is 
a pure state, its statistical properties are very different 
from the PSTD: for example, the TSRC does not exhibit 
squeezing for N > 12, as shown in the next section. 



III. STATISTICAL PROPERTIES OF TSRC 

A. Photon Number Distribution 

Since the expansion of TSRC is known in the number 
state \n), we have 



Pn =M 2 \C n \' 



(2) 
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Fig. 1 — 3 show the plots of the photon-number distribu- 
tion, P n versus n, for TSRC. Although the generation of 
the TSRC be virtually impossible for Hilbert space di- 
mensions 10 4 , due to nonidealities of beam splitters and 
detectors causing decoherence effects as well as due to 
low success probability (see section IV), it can be use- 
ful to study the behavior of TSRC for large TV, specially 
when we are interested in investigate the properties of 
sequence of numbers |14| . Fig— 3, which is a expanded 
view of the Fig. 2, shows the same behavior of that of 
the Fig.l. As expected, when N is large enough (Fig. 2 
), the difference from one peak to another becomes neg- 
ligible. For small values of n, running the RNG several 
times will give different sets of numbers with different 
P n s. Therefore, it is important to keep in mind that our 
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FIG. 1: The photon number distribution P n of TSRC versus 
n for the dimension of Hilbert space N = 10. This figure was 
obtained by a single run of the RNG. 



FIG. 3: Expanded view of the photon number distribution P n 
of the TSRC versus n in the Fig. 2. This figure was obtained 
by a single run of the RNG. 
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FIG. 2: The photon number distribution P n of TSRC versus 
n. As expected, when N is large the difference from one 
peak to another becomes negligible. The dimension of Hilbert 
space here is N — 10000. This figure was obtained by a single 
run of the RNG. 



B. Average number and variance 

The average number (n) as well as the variance in 
TSRC is obtained directly from 
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^P(n)n, 



(3) 



n=0 



and 



(An) 



(4) 



Fig. 4, 5 show the plots of (n) and (An) as a function 
of the Hilbert space dimension N. Note from Fig. 3 that 
(n) increases irregularly with N when a single realiza- 
tion is considered, but tends to increase linearly with N 
when several realizations are considered. As can be seen 
from Fig. 5, the increasing of (An) with N is irregular, 
irrespective of the number of realizations. 



C. Mandel parameter and second order correlation 
function 



Mandel Q parameter is defined as 



proposal to generate TSRC is based on a particular set 
of numbers obtained by a single run of the RNG. Thus, 
these numbers, having no algorithm producing them, a 
question emerging is: what type of properties will char- 
acterize a state having a random probability amplitude? 
These properties will be investigated in the next subsec- 
tions. 



(An 2 - (n)) 



(5) 



while the second order correlation function g^(0) is 
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FIG. 4: The average number < n > versus the dimension TV of 
the Hilbert space. Solid line was obtained from 30 realizations 
of the RNG. Dots refer to a single one. 
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FIG. 6: The Q parameter versus TV. For TV > 12, Q > 
and there is a linear dependence between Q and TV. Solid line 
was obtained from 30 realizations of the RNG. Dots refer to 
a single one. 
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FIG. 5: The variance (An) versus the dimension TV of the 
Hilbert space. 



and for Q < (Q > 0) the state is said to be sub- 
poissonian (super-poissonian) . Also, the Q parameter 
and the second order correlation function g^ are related 

by M 



Q 



,g (2) (0)-l (n) 



(7) 



It is well known that if g^(0) < 1 then the Glauber- 
Sudarshan P-function assumes negative values, thus dif- 
fering from usual probability distribution function. Be- 
side that, by Eq. (0) is readily seen that g^(0) < 1 
implies Q < 0. As for a coherent state Q = 0, a state is 
said to be a "classical" one if Q > 0. 

Fig. 6 and 7 show the plots of the Q parameter and the 
correlation function g^ 2 '(Q) versus TV. It is interesting to 
note that the TSRC is predominantly super-poissonian 
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FIG. 7: The g^ 2 \0) correlation function versus TV. Solid line 
was obtained from 30 realizations of the RNG. Dots refer to 
a single one. 



(Q > and g^(0) > 1) , then being a "classical" state 
in this aspect, for TV > 12. Nevertheless, for small values 
of TV (TV < 12), the Q parameter presents significant oc- 
currence of values less than 0, showing sub-poissonian 
statistics and thus being associated with a "quantum 
state". However, when TV becomes larger than 12, TSRC 
becomes super-poissonian. Also, note that for a single 
realization of the RNG these both functions present an 
oscillatory behavior, while, for several realizations, Fig. 6 
(solid line) shows an interesting linear dependence be- 
tween Q and TV. This will be further investigated in the 
following subsection. On the other hand, Fig. 7 indicates 
an asymptotic behavior of g^ 2 \0) when TV is increased. 
For TV = 10 4 , we find .g (2) (0) = 1.332585322. 
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FIG. 8: The AXi quadrature variance versus the dimension 
N of the Hilbert space. For N > 12, squeezing can be consid- 
ered a null event. Solid line was obtained from 30 realizations 
of the RNG. Dots refer to a single one. 



FIG. 9: The AX 2 quadrature variance versus the dimension 
N of the Hilbert space. For N > 12, squeezing can be consid- 
ered a null event. Solid line was obtained from 30 realizations 
of the RNG. Dots refer to a single one. 



D. Quadrature and variance 

Quadrature operators are denned as 

Xi = \ {a + o)); (8) 
X2 = ^(«-« t )- (9) 

where a (a') is the annihilation (creation) operator in 
Fock space. Quantum effects arise when the variance of 
one of the two quadrature attains a value AXi < 0.5, 
i = 1,2. Fig. 8 and 9 show the plots of quadrature vari- 
ance AXi = \J (X?) — (Xi) 2 versus N. Note the same 
oscillatory behavior when a single realization of RNG 
is considered. For small N (N < 12), TSRC exhibits 
squeezing with a significant frequency (running RNG 
several times). On the other hand, when N increases, 
squeezing becomes a rare event, and can be set a null 
event for N > 12, coinciding with the departure from 
sub-poissonian statistic discussed above. 
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FIG. 10: The S Shannon entropy versus the dimension N of 
the Hilbert space. Solid line was obtained from 30 realizations 
of the RNG. Dots refer to a single one. 



F. Husimi Q function 



E. Entropy 



Fig. 10 shows the plots of Shannon entropy S = 

N-l 

Pn^Pn versus N. It shows a steady increasing of 

n=0 

S. In fact, increasing N arbitrarily, we have found no 
bound for £, thus indicating no maximum for Shannon 
entropy. 



The Husimi Q-function for TSRC is given by 

Q\TSRC){P) = \\(m\ 2 ■ (10) 

Fig. 11 shows the Husimi Q-function for N = 15. It is to 
be noted that this function will change each time we run 
the RNG even for a fixed N. Thus, as remarked above, 
we refer to a particular fixed set of numbers assumed as 
having no law generating them. 
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FIG. 11: Husimi Q-function for the TSRC taking the dimen- 
sion of the Hilbert space TV = 15. 

IV. GENERATION OF TSRC 

Schemes for generation of an arbitrary state can be 
found in various contexts, as for example trapped ions 
[H and cavity QED HI HH- However, due to the se- 
vere limitation imposed by coherence loss and damping, 
we will employ the scheme introduced by Dakna et al. 
|17| in the realm of running wave field. For brevity, 
the present application only shows the relevant steps of 
Ref . |17| . where the reader will find more details. In this 
scheme, a desired state \^) composed of a finite number 
of Fock states \n) can be written as 

N N 

n— v ' n— 1 

r N 

= -^Y[D(f3 k )a + D((3 k )\0) 7 (11) 

* k— 1 

where D(f3 n ) stands for the displacement operator and 
the f3 n are the roots of the polynomial equation 

JY 

^ = 0. (12) 

n=0 

According to the experimental setup shown in the Fig.l 
of Ref. > we have (assuming 0-photon registered in all 
detectors) the outcome 

N 

|*) ~ J] D(a k+1 )a+T A D(a k ) |0) , (13) 
fe=i 

where T is the transmittance of the beam splitter and 
a k are experimental parameters. After some algebra, 



TABLE I: The roots (3* k = \/3 k \e~ iv ' 1 k of the characteristic 
polynomial and the displacement parameters a k = \ah\e.~ %lPo,h 
for TSRC. Here, N = 5, T = 0.878. The probability of suc- 
cessfully engineering TSRC is 0.2% 



N 


\Ph\ 




\0L k \ 


f ak 


1 


1.465 


3.141 


0.883 


0.297 


2 


1.080 


2.307 


0.647 


-2.316 


3 


1.080 


-2.307 


1.083 


1.570 


4 


2.190 


1.283 


2.465 


-2.005 


5 


2.190 


-1.283 


3.688 


1.570 


6 






2.190 


-1.283 



TABLE II: The roots /3£ = l/^e - ^ of the characteristic 
polynomial and the displacement parameters ol\ = \a k \e~ t{Pak 
for TSRC. Here, N = 13 and T = 0.951. The probability of 
successfully engineering TSRC is 10" 6 % 



N 


\h\ 




Kl 


<y3 Qfc 


1 


3.169 


3.141 


1.149 


0.102 


2 


2.651 


2.465 


1.089 


-2.156 


3 


2.651 


-2.465 


1.908 


1.570 


4 


1.529 


3.141 


1.055 


-1.885 


5 


0.701 


3.141 


0.526 


3.141 


6 


2.414 


1.731 


1.608 


-1.702 


7 


2.414 


-1.731 


3.353 


1.570 


8 


3.498 


1.420 


4.374 


-1.725 


9 


3.498 


-1.420 


5.381 


1.570 


10 


1.466 


0.921 


3.796 


-1.648 


11 


1.466 


-0.921 


2.009 


-1.648 


12 


2.772 


0.574 


2.747 


-2.065 


13 


2.772 


-0.574 


2.865 


1.570 


14 






2.772 


-0.574 



the Eq. ^H and Ea. (|13J) can be connected. In this 
way, one shows that they become identical when a\ = 
-EtiT-Wi and a k = T* N ~ k+l (f3 k -i - &) for 
fc = 2,3,4...7V. 

In the present case the coefficients C n are given by 
those of the TSRC. As stressed before, when N ^ 12, 
squeezing is lost and TSRC has its sub-poissonian statis- 
tics changed to a super-poissonian one. Thus, in the 
following, when considering values of N it can be inter- 
esting to take into account these two boundaries. Table 
I and II show the roots (3 k = \(3 k \ e~ l(pi3 k f the char- 
acteristic polynomial in Eq. l|12f) and the displacement 
parameters a k = \a k \ e~ lVa k for N = 5, and N = 13, 
respectively. The probability of successfully engineering 
TSRC, P\tsrc)i diminish with the number N, as is true 
for any state. 
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For N = 5, P\tsrc) = 0.2%. The beam-splitter trans- 
mittance which optimizes this probability is T — 0.878. 
For TV = 13, P\tsrc) — 10~ 6 % and the optimized beam- 
splitter transmittance is T = 0.951. 



V. FIDELITY OF GENERATION OF TSRC 

In the prior sections we have assumed ideal detectors 
and beam-splitters. Although very good beam-splitters 
are available by advanced technology (losses by absorp- 
tion lesser than 2%), the same is not true for photo- 
detectors in the optical domain. Thus, we now take into 
account the quantum efficiency r\ at the photodetectors. 
For this purpose, we use the Langevin operator technique 
as introduced in [2(j to obtain the fidelity to get the 
TSRC. 

Output operators accounting for the detection of a 
given field a reaching the detectors are given by |20| 



y/rja in + L a , 



(14) 



where 7/ stands for the efficiency of the detector and 
L a , acting on the environment states, is the noise (or 
Langevin) operator associated with losses into the detec- 
tors placed in the path of modes a = a,b. We assume 
that the detectors couple neither different modes a, b nor 
the Langevin operators L a ,Lf,. Thus, the following com- 
mutation relations are readily obtained from Ea. (|14(l : 



L L f 



L L f 



= l-r?, 
= 0, 



(15) 
(16) 



which give rise to following ground-state expectation val- 
ues for pairs of Langevin operators: 



(17) 
(18) 



ULH = l-r?, 



L L f 



0. 



These are useful relations specially for optical frequen- 
cies, when the state of the environment can be very well 
approximated by the vacuum state, even for room tem- 
perature. 

Let us now apply the scheme of the Ref.|l7J to the 
present case. For simplicity we will assume all detectors 
having high efficiency (77 > 0.9). This assumption al- 
lows us to simplify the resulting expression by neglecting 
terms of order higher than (1 — rf) 2 . When we do that, 
instead of the ideal \TSRC), we find the mixed state 
\^fe) describing the field plus environment, the latter 
being due to losses coming from the nonunit efficiency 



detectors. The result is, 

\V FE ) ~ R N D(a N+1 )a + T A D(a N )a+T A 
x D{a N ^ 1 )...a + T il D{ai) |0)L+ 
+ R N - 1 D{a N+1 )a+T h D{a N )a + T fl 
x D(a N ^ 1 )..X+T A D(a 1 )\0) 



R 



N- 



1 I?(aAr +1 )a + T fi ^(a w )L+_ 1 T fl 

x D(a N ^ 1 )...a + T h D(a 1 ) |0) 

+ R N - 1 D{a N+1 )L+T fl D(a N )a + T fl 

x D(a tf -i)...d+T fi D(a 1 ) |0) , (19) 

where, for brevity, we have omitted the kets correspond- 
ing to the environment. Here R is the reflectance of the 
beam splitter, = 1 is the identity operator, and L^, 
k = 1, 2..N stands for losses in the first, second . . . N-t/i 
detector. Although the L^/s commute with any system 
operator, we have maintained the order above to keep 
clear the set of possibilities for photo absorption: the 
first term, which includes Lq = 1, indicates the proba- 
bility for nonabsorption; the second term, which include 
Lf , indicates the probability for absorption in the first 
detector; and so on. Note that in case of absorption at 
the k-th detector, the annihilation operator a is replaced 
by the L^" creation Langevin operator. Other possibili- 
ties such as absorption in more than one detector lead 
to a probability of order lesser than (1 — 77) 2 , which are 
neglected. 

Next, we have to compute the fidelity jSJ, F = 
|| |\Pf.e)|| , where I*) is the ideal state given by 
Ea. (|13|) . here corresponding to our CCS characterized 
by the parameters shown in Table 1, and \^fe) is the 
(mixed) state given in the Ea. l(T§|) . Assuming 77 = 0.950 
and 0.900 we find, for N = 5 (13), F = 0.997 (0.994) 
and 0.996 (0.991), respectively. These high values of fi- 
delities show that efficiencies around 0.9 lead to states 
whose degradation due to losses is not so dramatic. 



VI. COMMENTS AND CONCLUSION 

In this paper we introduce a new state of the quantized 
electromagnetic field, the Truncated State with Random 
Coefficients (TSRC). To characterize the TSRC, we have 
studied several of its statistical properties as well as the 
behaviour of these statistical properties when the dimen- 
sion N of Hilbert space is increased. By studying the de- 
pendence of the Q-parameter and of the second order cor- 
relation function with increasing N, we found that there 
exists a transition from sub-poissonian statistics to super- 
poissonian statistics when N is relatively small (N ~ 12). 
Interesting, we have also found that for N > 12 squeez- 
ing is rarely observed. Other statistical properties such 
as Shannon entropy, quadrature operator and their vari- 
ances were investigated. We think TSRC can be used to 
characterize sequence of numbers generated determinis- 
tically. It would be interesting, for example, to compare 
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TSRC properties with the properties of states whose co- CAPES (WBC) and the CNPq (NGA), Brazilian agen- 
cfhcients stem from a sequence generated by iteration of cies, for the partial supports of this work, 
the logistic equation. 
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